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Abstract 

Specific heat measurements down to 0.5 K have been performed on a single 
crystal sample of a spin-ladder like compound (^(CsH^^^CU under mag- 
netic fields up to 12 T. The temperature dependence of the observed data in 
a magnetic field below 6 T is well reproduced by numerical results calculated 
for the 5=1/2 two-leg ladder with </ rU ng/«/ieg=5. In the gapless region above 
7 T (-Hci)j the agreement between experiment and calculation is good above 
about 2 K and a sharp and a round peak were observed below 2 K in a mag- 
netic field around 10 T, but the numerical data show only a round peak, the 
magnitude of which is smaller than that of the observed one. The origin of 
the sharp peak and the difference between the experimental and numerical 
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round peak are discussed. 
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I. INTRODUCTION 



Recently, there has been a considerable interest in quantum spin systems with a spin gap 
above the singlet ground state. One of the examples studied extensively are one dimensional 
Heisenberg antiferromagnets (lDHAFs) with integer spin values, especially spin(5) one, 
which is associated with Haldane's prediction M. Now, these studies extend those of the 
5=1 antiferromagnetic bond alternating chains Another case is 5=1/2 two leg spin- 

ladder which is investigated as an intermediate system between one- and two-dimensional 
systems [UQ. Spin- ladder systems have been studied in relation to the Haldane problem 
on one hand @J^] and high T c superconductivity on the other hand PH. Most of the 
spin-ladder systems investigated so far are copper oxides pJlQfl. One of remarkable things 
is that a lightly hole doped two leg ladder exhibits superconductivity under high pressure 
PJTT|] as expected theoretically [[§]. On the other hand, spin-ladder like copper complexes 
Cu2(C5H 12 N 2 )2Cl2( abbreviated as CHpC) [|T2| — pL9|| and KCuCl 3 and its family compounds 
22) have been also studied extensively. These compounds except for NH4CUCI3 f2~I[ 



have the singlet ground state and become gapless at a certain magnetic field (H&). The 
features of this gapless region at low temperatures have attracted much interest because the 



field induced long range ordering (LRO) in TICUCI3 |22j has been interpreted as a Bose- 
Einstein condensation of magnons |2"3|] . This matter was originally argued by Affleck [23 
that the ground state above H c \ may be regarded as a Bose condensate of the low energy 
boson. Besides, the dimensionality is expected to appear in the power-law dependences of 
thermodynamic quantities on temperature, when approaching the quantum critical point by 



application of a magnetic field | 25fl . Furthermore, in the quantum critical region above the 
LRO temperature, we can investigate the feature of Tomonaga-Luttinger(TL) liquid for a 
quasi one- dimensional antiferromagnet • 



The compound CHpC has been studied extensively for several years by various exper- 
iments such as magnetic susceptibility Hj- ljj] , magnetization [12|-p|, NMR |15,|16|, ESR 
I1] , |PT|| , specific heat |fl2|JT8| , p!9[] and neutron scattering measurements jnj. From these ex- 
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periments, CHpC has the singlet ground state with an excitation gap of about 10 K and 
the gap collapses at about 7 T (H c i) and the saturation field is about 13 T (H c2 ). This 
compound has an advantage of study on the gapless spin liquid phase under magnetic fields, 
because H c i and H c2 are easily accessible fields with a conventional superconducting mag- 
net. Specific heat measurements in a magnetic field were done up to 9 T by Hammer et 
al. Jl^l and up to 8.25 T by Calemczuk et al. |18]. In the present paper, we extend the 



specific heat measurements under magnetic fields up to 12 T beyond the symmetry field 
H 8yril (—{H c i+H C 2) /2) and compare the experimental results with those of numerical and 
analytical calculations. 

The format used in this paper is as follows: In Sec. II, experimental and theoretical 
details are described. Experimental results of magnetic susceptibility and specific heat 
under magnetic fields(if ) are reported in comparison with some numerical calculations in 
Sec. III. A sharp and a round peak observed at low temperatures in the gapless region are 
discussed in Sec. IV. 

II. EXPERIMENTAL AND NUMERICAL DETAILS 



Powder samples of CHpC were synthesized according to the method reported in Ref.[26]. 
Equimolar amounts of 1,4-diazacycloheptane (C5H 12 N 2 ) and CuCl2-2H 2 were dissolved in 
warm methanol (60° C) for 1 hour and left at room temperature for two days. Single crystals 
of CHpC were obtained by the slow evaporation method from a methanol solution of powder 
samples of CHpC. We obtained samples with 2x2x1 mm 3 in typical size. 

CHpC crystallizes in the monoclinic system and belongs to the P1\jc space group ||28|| . 



The lattice constants and (3 angle at room temperature are a=13.406(3)A, 6=11.454(2)A, 
c=12.605(3)Aand /5=115.01(2)°. Cu dimeric units are linked to the neighboring ones via 
hydrogen bonds along the [101] direction. The upper panel of Fig. 1 shows a schematic view 
of the chain like structure along [101] of CHpC. Broken lines in the upper panel show the 
hydrogen bonds and the ellipsoids around the copper atoms represent the 3d hole orbitals. 
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In the lower panel of Fig. 1, possible pathways of the exchange interaction between the 
neighboring Cu 2+ spins are depicted. The exchange interaction on the broken pathway must 
be much weaker than those on the other pathways because of the configuration of 3d copper 
and 3p chlorine orbitals, thus considering this system as the S—l/2 two leg spin-ladder. 

Magnetic susceptibilities (M/H) were measured with a SQUID magnetometer (Quantum 
Design's MPMS2) in RIKEN. Specific heat measurements down to 0.5 K under magnetic 
fields up to 12 T were performed with Mag Lab^ c micro calorimeter (Oxford Instruments) 
installed at the same place. The relaxation method was employed so that we used only one 
single crystal with 4.3 mg in weight. Numerical calculations were done by the temperature- 
dependent density- matrix- renormalization-group (DMRG) method for the S—l/2 two- leg 
Heisenberg spin-ladder with a spin Hamiltonian written as 

N/2 2 N/2 

H. = Jrung X ■ S 2 J + Jleg ^ X ' Sjj+i - Q^H ^ S-j, (1) 

j=l i=l j=l i,j 

where J rung and J\ eg are the exchange constants along the rung and the leg, respectively 
and Sij an S—l/2 spin on the i-th leg, the j-th rung, g the g-value of copper, /ie the Bohr 
magneton and H the external magnetic field. As we use a transfer matrix approach, it is 
free of finite size effects. 

III. EXPERIMENTAL RESULTS AND COMPARISON WITH CALCULATIONS 

A. Magnetic susceptibility 

Solid circles of Fig. 2 shows the temperature dependence of magnetic susceptibility of 
a single crystal sample of CHpC along the chain direction. A hump which is typical of 
low dimensional antiferromagnets is observed around 10 K and the susceptibility steeply 
decreases with further decrease of temperature toward zero Kelvin. No increase of the 
susceptibility due to magnetic impurity or crystal defect is observed in this sample. The 
solid line in this figure is a fit of the calculated susceptibilities to the experimental ones with 
the fitting parameters of g—2.1, J mng =13.1 K and Ji eg =2.62 K ( J run g/ •7ieg=5). The agreement 
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between numerical and experimental susceptibilities is excellent. Evaluated values and the 
ratio of the exchange constants are close to those estimated by other groups P2,I3|. In the 
following comparison of specific heat, we use the same fitting parameters. 



B. Specific heat at H=0 T 

In Fig. 3, filled circles represent raw specific heat data. Usually the lattice part of the 
specific heat is evaluated from the specific heat of a nonmagnetic isomorphous compound. 
But no isomorphous compound with nonmagnetic atom like Zn exists. In order to get the 
magnetic contribution of the specific heat, we subtract the lattice part (Ci att i ce ) of the specific 
heat from the raw data in the standard way, assuming Ci at tice~^ 3 at low temperatures. The 
broken line in Fig. 3 is the lattice part of the specific heat with a coefficient of 0.004 (J/K 
Cu-mol) which is not far from those evaluated in other copper complexes (e.g. 0.0048 for 



(CH3)2CHNH 3 CuCl3 |29|) . Open circles show the subtracted results, namely, magnetic 
part of the specific heat and we see a round peak due to the short range ordering around 
5 K. The solid line represents the specific heat at H—Q T calculated with the same fitting 
parameters as the susceptibility fitting. The calculated specific heat satisfactorily agrees 
with the observed one. At the temperature below the round peak, an exponential decay due 
to a spin gap is observed. We fit the specific heat data at low temperatures to an expression 
of C mag ~T~ 3 / 2 exp(-A/T) to evaluate the energy gap. This expression is deduced in the low 
temperature limit for the strong coupling case ( J r ung/^ieg=5) approximating the dispersion 
J rung + JiegCosg where q is the component of wave vector transfer along the chain [fT^| . The 
fitting result is drawn as the solid line in Fig. 4 with the fitting parameter A=10.9 K, which 
is very close to the value of J rung - Ji eg (=10.48 K). 

C. Specific heat at H ^0 T 

Filled symbols in Figs. 5(a), 5(b) and 5(c) show the magnetic specific heat data for 
magnetic fields below H c \, between H c \ and H sym and above H sym , respectively. In Figs. 5(a), 
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5(b) and 5(c), a hump is observed above 5 K at each designated magnetic field. Moreover, 
we observed two peaks below 2 K near H sym in Figs. 5(b) and 5(c). Details of this low 
temperature part are depicted in the upper panel of Fig.6. We obviously see two peaks, a 
sharp and a round peak at the magnetic field above 8.5 T. In this figure, we plot the specific 
heat data under 10. 5T for simplicity. Figure 7 shows the magnetic field dependence of the 
specific heat. Two peaks are observed below 0.82 K (only one peak at 0.61 K because of our 
instrumental limitation), whereas no peak is observed at 0.87 K. We plot these two peaks 
at low temperatures in Figs. 5-7 in the plane of H vs. T of Fig. 8. Plotted points are almost 
symmetric at if sym (~10 T) for both peaks. 

Next, we compare the experimental data with the numerical ones. Solid, broken, dotted 
lines in each figure of Fig. 5 represent the numerical specific heat data for the designated 
magnetic fields. In Fig. 5 (a), the agreement between experimental and numerical specific 
heat is excellent over the entire temperature range up to 9 K. In Figs. 5(b) and 5(c), the 
agreement between experiment and calculation is good above 2 K, while the large deviation 
is seen below about 2 K. In the lower panel of Fig. 6, calculated specific heat data are 
displayed for the magnetic fields corresponding to those in the upper panel. The tendency 
of the calculated round peak to shift with change of fields is similar to the observed one, 
but the magnitude of the peak is much smaller than that of the observed peak. We show 
in Fig. 9 the magnitude of the round peak as a function of magnetic field. The field at the 
calculated maximum peak slightly shifts to the lower side compared to the experimental 
one. The magnitude of the calculated peak is almost a half of the observed one. 

Entropy of this sample is calculated from the specific heat data shown in the upper 
panel of Fig. 10 for the designated magnetic field. Correspondingly, we show the calculated 
entropy for the corresponding magnetic field in the lower panel of Fig. 10. Similar tendency 
of the entropy is observed in both figures. With increasing the magnetic fields, the entropy 
becomes higher at low temperatures. This behavior is opposite to that in a paramagnet. 
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IV. DISCUSSION AND CONCLUSIONS 



First let us discuss the origin of the sharp peak observed in the experiment. Usually, 
the sharp peak is thought to be caused by magnetic long range ordering (LRO) [^7j, but 
it is controversial for the case of a spin ladder; for example, it predicts a field dependence 
of the transition temperature (camel-type structure) contrary to experimental observations 
(dromendary-type structure). Recently, Nagaosa and Murakami [RO] argued that a lattice 
instability in the spin ladder is expected to occur above H c \. From their study, the lattice 
distortion occurs in the spin ladder at an incommensurate wave vector corresponding to the 
magnetization. 

For a comparison of these scenarios with our experiment, we note that the one of the 
experimental consequences of Nagaosa's and Murakami's scenario is that a gap should appear 
below the transition temperature (T c ) and the specific heat should behave as C ~exp(- 
jv /T), where jv is the gap of the order of T c . In the upper panel of Fig. 11, we fit the 
experimental data at 10 T below the temperature at the sharp peak to this equation and 
obtained fairly good agreement between them with 7t>o=2.58 K, which is however quite far 
from the temperature T c of about 0.8 K at 10 T. If the sharp peak is on the other hand caused 
by antiferromagnetic long range order, the magnetic part of the specific heat should go as 
T 3 due to the antiferrmagnetic magnons. We show the fitting result using this equation in 
the lower part of Fig. 11. The agreement between experiment and calculation is comparably 
good to the previous fitting, when we add a constant. This negative contribution is however 
unexplained. Both interpretations are therefore not really conclusive, and discrepancies 
remain between the possible sources of the sharp peak and the specific heat data, calling for 
the proposal of a new transition scenario. 

Theoretically we can discuss the round peak based on a fermion representation. The 
round peak can be attributed to the low energy excitation structure in a gapless spin liquid 
(TL liquid). As pointed out in the related context of mixed spin 1 and \ by Kolezhuk et.al. 
PHI , the round peak arises from a spinon band in the spin liquid system. They reduce the 
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Hilbert space of the problem, keeping only the most important states per elementary cell. In 
the strong coupling case of the spin ladder, this amounts to keeping for each rung only the 
singlet and the lowest-energy triplet component around H c i, mapping them to an effective 



spin 1/2 chain |2E|J27]. Then the effective Hamiltonian is given by 



N N 

W = J* E («? • s t+i + s " ■ s li + V2a? " - {g^H - J rung - l/2J leg ) £ «?, (2) 

i=l i=l 

where sf, and sf are the x, y and 2 component of the effective 1/2 spin. When, for a 
qualitative consideration, the sf ■ s z i+1 term is neglected, the remaining free fermion hamil- 



tonian can be easily treated. Then a splitting of the peak is expected to occur |31[ due to 
contributions from two different spinon bands (particle and hole band) for magnetic fields 
off the symmetry field H sym . In experiment and in numerical studies on the full original 
Hamiltonian of the ladder, we however observe only one peak, which is probably due to the 
large deviation of the effective Hamiltonian from the exactly solvable xy-limit, making the 
analytical prediction of a peak split less stringent. 

Finally, the difference of the magnitude between the experimental and the numerical 
round peak has to be discussed. When calculating the specific heat, we neglected the 
diagonal exchange interaction denoted by j in Fig. 1. This interaction can possibly affect 
the nature of the specific heat at low temperatures, but it should not allow for an effect 
of more than 50 percent as observed. There must be an additional source of entropy and 
disorder, probably related to the physics driving the three dimensional transition, which at 
present can only be speculated about. 

In conclusion, we have made specific heat meaurements on a single crystal sample of the 
S—l/2 Heisenberg spin-ladder compound Cu2(C5Hi 2 N 2 )2Cl2 under magnetic fields. Exper- 
imental data are compared to the numerical ones calculated with a temperature dependent 
DMRG method for the 5=1/2 spin-ladder and are well reproduced when H < if cl , but when 
H > H c i and for temperatures below about 2 K the magnitude of the observed round peak 
differs by a factor of about 2. The numerical calculations give quasiexact results for the ID 
Hamiltonian considered in the whole temperature and magnetic field range. The remaining 
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discrepancies have therefore to be due to effects beyond that Hamiltonian, probably due to 
low-temperature 3D coupling and/or further degrees of freedom (such as phonons), as also 
indicated by the 3D transition. 

We have observed two peaks, a sharp and a round peak at low temperatures below 2 K 
around H sym . The origin of the former peak is still controversial and two possibilities are 
discussed: the antiferromagnetic LRO |27fl and the lattice instability above H c \ proposed by 



Nagaosa and Murakami [p0| . The latter peak probably comes from a spinon band in the 
effective Tomonaga-Luttinger liquid. A simple analytical fermion representation predicts 
moreover a weak splitting of the round peak into two round peaks for magnetic fields off 
the symmetry field H Bym , but experimental and numerical results show consistently only one 
peak. This shows that on the analytical level a more elaborate analysis of the Hamiltonian 
is needed to describe the round peak in the ID picture, while qualitatively the emergence 
of a new low-temperature structure in the specific heat is well captured. 
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FIGURES 

FIG. 1. A schematic view of the chain structure of (^(CsH^^^CU and the exchange path- 
ways. Broken lines represent hydrogen bonds and ellipsoides show 3d hole orbitals of copper. 

FIG. 2. Temperature dependence of magnetic susceptibilities of a single crystal of 
Cu2(C5Hi2N2)2Cl4 along the chain direction ([101]). The solid line represents calculated sus- 
ceptibilities for the 5=1/2 Heisenberg spin-ladder with the fitting parameters shown in the panel. 

FIG. 3. Specific heat as a function of temperature at H=0 T. Filled circles are raw specific heat 
data and a broken line represents the estimated lattice part of the specific heat. Open circles show 
the magnetic specific heat. The solid line represents the specific heat calculated for the 5=1/2 
spin ladder with the same fitting parameters as the susceptibility. 

FIG. 4. Specific heat at H=0 T in the plane of ln[C mag T 3 / 2 ] vs. 1/T. The solid line is the 
result of a fit to the equation ~T 3 / 2 exp(-A/T) with A=10.9 K. 

FIG. 5. Temperature dependence of the specific heat under the designated magnetic fields 
along [101] direction. The solid, broken and dotted lines are calculated specific heat data for the 
corresponding magnetic fields. 

FIG. 6. Details of the specific heat as a function of temperature at low temperatures below 1.5 
K above 

FIG. 7. Magnetic field dependence of the specific heat along the [101] direction for the desig- 
nated temperatures. Solid lines are guides for the eye. 

FIG. 8. Plot of the peaks of the specific heat at low temperatures above 7 T in the magnetic 
field vs. temperature plane. Filled circles and open squares correspond to the sharp peaks and the 
round ones, respectively. Solid lines are guides for the eye. 
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FIG. 9. Magnitude of the round peak as a function of magnetic field. Filled circles and open 
squares represent the magnetic specific heat data for experiment and calculation, respectively. Solid 
lines are guides for the eye. 

FIG. 10. (a) Entropy as a function of temperature calculated from the observed specific heat 
data for designated magnetic fields. The entropy at the lowest temperature of about 0.5 K is 
evaluated as half of C mag at the lowest temperature, (b) Calculated entropy as a function of 
temperature corresponding to the experiment. Solid lines are guides for the eye. 

FIG. 11. (a) Specific heat data at 10 T in the plane of -lnC mag versus 1/T. The solid line 
represents a fit of the experimental data to an equation ~exp(-7Wo/T). (b) Specific heat data at 
10 T in the plane of C mag versus T 3 . The solid line shows a fit of the experimental data to an 
equation a+bT 3 where a and b are fitting parameters. 
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Fig. 1 M. Hagiwara et al. 
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